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Pumping water is one of the most popular technologies of solar energy for irri- 
gation or drinking water supply. Its performance depends on the characteristics 


of the site (sunlight, ambient temperature, geographical obstacles), on the per- 
formance of the modules, as well as, on the characteristics of other equipment 
(converter and pump). The optimal sizing of a photovoltaic water pumping sys- 
tem makes it possible to guarantee the satisfaction of the water demand through- 
out the irrigation period and to model the electrical energy needs to supply the 
pump to irrigate the crops and water livestock. This paper proposed a novel al- 
gorithm for dimensioning the elements of an autonomous photovoltaic system 
equipped with irrigation water storage. The proposed algorithm determines the 
optimal surface area of the photovoltaic modules and the electrical power nec- 
essary to satisfy the water requirement for irrigation in the observed time. The 
results obtained from the proposed algorithm for optimal sizing are compared to 
calculation scenarios to show better results. A case study from the Fez-Meknes 
region in Morocco has been selected for applying the optimal sizing algorithm. 
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1. INTRODUCTION 

Photovoltaic solar pumping over the sun is an efficient solution, in terms of energy, for the supply of 
water in the absence of the electricity grid. This type of system is well suited to applications where a constant 
volume of water must be pumped daily from a water source to the application or storage of water, as it’s well 
simply because it consists of a direct coupling between the photovoltaic generator and the pump, and it is 
intended for pumping time throughout the day. In this case, the system can store energy or water to store 
the excess energy generated on days when solar energy is above average, either by storing extra electricity 
generated in batteries [1], or by pumping extra water and storing it. 

One of the major issues facing solar irrigation systems is the seasonal demand pattern for crops and 
the diversity of solar power generation. To cope with this problem, several photovoltaic systems use batteries to 
accumulate the electrical energy produced by the photovoltaic system to match energy production and demand 
[1]. However, the batteries have several disadvantages because they are heavy, expensive, fragile, and require 
maintenance of the battery. Therefore, batteries are also introducing some degree of efficiency loss from about 
20% to 30% of energy production [2], it causes increased system investment costs, space to store batteries, and 
environmental concerns related to battery disposals [3]. 
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Several approaches have been proposed for a self-contained, direct-pumped photovoltaic irrigation 
system that pumps water directly to the irrigation distribution network instead of storing it in a reservoir. 
Moreover, to match energy production with irrigation demand [3], [4]. Adapting the production of photovoltaic 
energy to the demand for crop irrigation and also the livestock’s demand for storing pumped water in reservoir 
storage remains the most relevant solution in terms of the availability of water in the worst monthly average of 
solar radiation. 

The simplicity of the solar photovoltaic pumping system faces the problem of load adaptation since 
the direct coupling doesn’t allow the generator to deliver its maximum power for a full day. It is, therefore, 
necessary to operate these generators at their optimum power. The approach to the problem is to insert between 
the generator and the receiver a static converter (chopper or inverter) to perform the optimal transfer of power. 
Furthermore, direct-pumping stand-alone photovoltaic irrigation systems that pump water directly to the asking 
location can cause problems where the water demand will not be sufficient in the critical time. This is why we 
have proposed a system to combine direct irrigation and irrigation through the storage tank. This solution meets 
the need at any given moment. 

The main objective of this work is to develop a novel algorithm to provide optimal dimensions for 
the components of photovoltaic (PV) pumping systems to ensure that the water demand is met throughout the 
irrigation duration. This algorithm was applied to stand-alone PV pumping systems to irrigate tomato farmland 
located in the Fez-Meknes region of Morocco during the most critical month of the culture cycle vegetative, 
which runs from April to August. The system consists of PV modules, a DC/DC converter, livestock watering, 
and a water tank. Figure 1 shows the principal components of the standalone photovoltaic system used to pump 
water in this application. 
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Figure 1. Synoptic scheme of stand-alone PV pumping system 


This paper presents three sections. The first section (section 2) is the modeling of the parts of electrical 
and hydraulic energies: solar radiation and global irradiation calculation model, optimal position of the panel, 
photovoltaic array model, and motor-pump unit model. The second (section 3) is about the water requirement 
for irrigation, livestock, and irrigation water storage sizing. The third section (section 4) presents the methodol- 
ogy of the proposed algorithm. Then, a case study to test the performance of the sizing algorithm and compared 
results which are described and explained in section 5. The conclusion is given in section 6. 


2. SYSTEM MODELLING 
To dimension the elements of the system, it is necessary to model the components of the installation. 
This modelling includes both electric and hydraulic parts which are presented in this section. 


2.1. Photovoltaic modules 

Global solar radiation is a function of the composition and thickness of the atmosphere traversed by 
the light rays during the day. The total daily solar radiation G.,(a, 8) in a tilted panel is evaluated by changing 
the hour angle corresponding to the length of the day. It is expressed by (1) [5]. 
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(1 + cos) 


Gs(a, 8) = Raln + Dp (1 S a 


+ Alb Gh (1) 
Where J}, is the direct normal irradiance, Dp is the diffuse horizontal irradiance, Ghis the global horizontal 
irradiance, Alb is the albedo of the soil, 3 is the module tilt angle and Rg is the ratio of direct radiation on the 
inclined panel and direct radiation on the horizontal panel. 

To model a photovoltaic pumping system, it’s first necessary to model the photovoltaic cells. There 
are several mathematical models [5], [6] like one-diode and two-diode models, which describe the photovoltaic 
current as a function of the photovoltaic voltage, the multi-parameter I-V which is a non-linear equation. These 
photovoltaic models cell take into account the influence of the daily variation of global radiation and the junc- 
tion temperature. The difference between these models is the calculation method, the number of intermediate 
parameters to calculate the I-V characteristic, and the accuracy of the results. 

The model used in this paper is a single model (Figure 2), which is an analytical model that has given 
a great accuracy and can be suitably used for power system planning purposes [7], [8]. The I-V characteristic 
of the PV module is given by [5]-[9]. 


(2) 


Ve + IRs Ve + Rsle 
I, am lon = I, (exx Es ) 1) ( bs ) 


V: Rp 


Where Rp is the parallel resistance of the PV module and Rs is the serial resistance of the PV module. The 
panel yield model is expressed by [5]-[10]. 


Mpv = Nr [L — bpv (Te — Tres )] (3) 


Where ny is the panel yield at the reference temperature, py is the temperature coefficient for the panel yield, 
Te is the cell temperature, and T,., ¢ is the reference temperature. The cell temperature is varying in the function 
of solar radiation and the ambient temperature. It can be calculated as [11]. 


NOCT — 20 
800 
Where T, is an ambient temperature, NOCT is the Nominal Operating Cell Temperature and G's (a, 8) is the 
solar radiation on the tilted panel. 7,, 8p, et NOCT depend of the type of photovoltaic array considered. 

Eventually, the PV power can be assessed by [12]. 


Te = Ta + ( )+ Gs(a, B) (4) 


Poy = S- Gs(a, B) > Nye (5) 


Where S is the panel surface. 


Figure 2. Photovoltaic cell model 


The sun turns from the east, where it rises, to the west, where it sets every day. This varies during the 
seasons. It varies all the more as the latitude of the place of observation is high. During the summer season, the 
sun rises to the south-east of the sky and lies to the south-west. On the other hand, during the winter season, 
the sun rises approximately northeast of the sky and sets to the north-west. So, in order to maximize the power 
output of a photovoltaic system, it is necessary to orient the modules optimally in order to capture a maximum 
of solar radiation. Indeed, the direction of the sun is important. 
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There are many methods of sun-tracking that have been proposed to facilitate the tracking of the sun’s 
position with a high degree of accuracy [13]-[16]. The panels can absorb the maximum energy in the noontime, 
where the rise of the sun is at zenith and the values of the solar radiation are maximum. Furthermore, for the 
trackers and in order to meet the energy needs and optimize the overall efficiency of the system, the sensor field 
must keel its surface perpendicular to the rays emitted by the sun. 

To calculate the optimal position of the panel described by the two angles (inclination 8 and orientation 
a), we must look for these two angles from the radiation incidence angle (6), for which the angle of incidence 
y will be zero (8 = 90° — a) and the module azimuth angle equal to the sun’s azimuth angle (a = a). 


cos(y) = cos(8)sin(h) + sin(B)cos(h)cos(a — a) (6) 


Where h is the sun elevation angle, a is the surface azimuth angle and a is the solar azimuth angle. 


y = 0° > cos(y) = 1 
B=90° -a ) 


a=a 


In the practical case, it’s very difficult to implement these methods because of the number of panels 
and the exorbitant cost generated depending on the mobile support. In addition, they will have to be maintained 
regularly and may be more likely to be damaged by high winds; so fixed systems are generally cheap and 
require the least maintenance. In these methods, the solar panels are installed at an inclination 8 and a fixed 
orientation a. The latter is chosen towards the south or the north, according to the geographical situation of 
the considered site compared to the equator. For the case of Morocco, it is located in the northern hemisphere, 
so the orientation of the panels is towards the south. For inclination, the modules must be tilted at the latitude 
+15°(+ for the summer, - for the winter) [17] of the site to produce the largest annual solar radiation and hence 
the highest annual output power at that location. 

Photovoltaic modules could be oriented to the east or west to maximize energy production in the 
morning or afternoon as needed. It is also possible to tilt the photovoltaic modules at an angle greater than the 
latitude angle if the need for water pumping is higher in winter or throughout the year. On the other hand, they 
could be inclined at an angle less than the latitude angle if the pumping requirements are higher in summer. 


2.2. Pump modelling 

Each pump features a flow rate and total head. The typical head consists of a static component and 
a dynamic component [18]. The static head is the vertical distance from which the water must be pumped. It 
is the height difference in feet between the pumping level in the well and the pressure tank (static head = H, 
+ H.,). The total dynamic head (T DH) is the sum of the static head, the draw-down distance (the water level 
drops), and the distance equivalent to the friction losses in the pipe; it can be expressed by [19]. 


TDH = H.+ H, +H, + Fi (8) 


Where He is the vertical rise, (H, + H,) is the pumping level and fF; is the friction loss. 

A larger pipe size must be uses in photovoltaic pumping systems than would be used with conventional 
water systems to maintain friction losses as low as possible [20]. It is important to match the discharge capacity 
of the well and the pumped water, i.e. the power of the PV pumping system must be synchronized in this sense. 
As a result, it is necessary to introduce a limitation on the daily pumping of water from the well in the period 7, 
with the following equation [4]-[12]. 

Q(t) < Qmacti (9) 


Where i assumes the values 1 to Nr (Nr is the total number of time stages, decades), Q(i) is the mean daily 
pumped water (m°/day), Qmazx is the maximum discharge capacity of borehole (m/h), t; is the mean daily 
insolation (h). 


3. WATER REQUIREMENT 
Determining the total amount of water needed per day should be the first step in sizing the system. In 
our case, we have two parts to determine its need for water: irrigation and livestock watering. 
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3.1. Crop water requirement 
The crop water requirements correspond to how much water to bring to the plant. Water requirements 
vary mainly depending on the type of plant, soil, climate, and weather conditions. The balance between the 
inputs (rain, watering, groundwater, and the residual moisture of the soil) and the losses which are modelled 
in the form of evapotranspiration is made to know the quantity of water to be brought. It can be represented 
by [21]. 
V, = ET. — P, — ASF — AW (10) 


Where V,, is the net water requirements, ÆT, is the culture evapotranspiration, P. is the effective precipitation, 
ASF is the water which arrives at the root zone by capillary ascent from the phreatic level, and AW is the 
variation of residual soil moisture. 

According to [22], evapotranspiration is a parameter that depends on the climate, the type of plant, 
and its state of development. It can be expressed as (11). 


ET, = K.ETo (11) 


Where Ke is the crop coefficient, which depends on the characteristics of the crop itself (growth period, degree 
of water supply of crops during the vegetation period, and climate conditions of the area), and FT is the 
amount of water that plants may use from the soil and the surface of plant organs under optimum cultivation 
conditions [21]. It is essential to estimate the irrigation losses to obtain the total water requirement, which 
represents: i) evaporation losses at ground level; ii) losses by runoff or deep percolation; iii) soil washing, 
when more water is needed to keep the salts in the soil away from the roots; and iv) inhomogeneous water 
distribution, which means bringing more water to the most favorable areas to ensure that the most unfavorable 
areas receive a sufficient amount. For this, total water requirements during the time stage evaluated by [21]. 


_ ET. — Pz 


Vri 
We 


(12) 


Where We is the watering efficiency. 


3.2. Livestock requirement 

The water requirements for livestock depend on the type and quantity of animals. These needs vary 
with air temperature, age, herd size, gestation, breastfeeding, animal weight, type of diet, and size of livestock. 
It should be taken all of these variants into account when determining the minimum volume of water that the 
solar pumping system must provide each day. Water requirements correspond closely to the amount of food 
or feed ingested; as intake increases, the level of water requirements increases [18]. Hence, the daily water 
requirement per day (Gw) can be calculated as [23]-[25]. 


Gu/Item = Qz * R,/Item (13) 


Where Qz is the quantity and Rọ is the required gallons of water per day which is 3.78541 liter of water per 
day. 


3.3. Irrigation water storage 

Pumping over the sun makes it possible to have a photovoltaic system that is simpler, more reliable, 
and less expensive than a system with a battery. A water storage tank is an essential component of a solar- 
powered water pumping system and is economical. The storage is done hydraulically, the water being pumped, 
when there is enough sunlight, and the temperature of the solar cells, in a tank above the ground. It is then 
distributed by gravity as needed (livestock and irrigation). The reservoir can often be built locally, and the 
storage capacity can vary from one to several days. This tank doesn’t require complex maintenance and is easy 
to repair locally. For stand-alone systems, the storage tank must have sufficient volume to accommodate peak 
demands, to compensate for cloudy days and nights. For this reason, the required volume of the storage tank is 
determined by (14) [5]. 


Vstorage = Nstorage(Vi = Vri) (14) 


Where 1)storage is the water losses in the storage tank, V; is the water volume leaked or in excess, and V,; is the 
total water requirements for irrigation. 


A novel algorithm for optimal sizing of stand-alone photovoltaic ... (Asmae Hafian) 


1838 m) ISSN: 2088-8694 


4. PROPOSED ALGORITHM 

The objective is to provide optimum component dimensions that can provide autonomy for the pump 
installation and pump a volume of water needed to irrigate crops and livestock. This algorithm depends on 
panel and site characteristics, including solar radiation, ambient temperature, and crop characteristics. The 
algorithm flowchart is presented in Figure 3. The algorithm is performed by Algorithm 1. 


Start 


| 


Panel, Site and Crop characteristics 


| 


i. Estimation of the diffused and direct v. Evaluation of the total water 
radiation from the site characteristics; requirements for irrigation Vri using 
Eq.12; 


ii. Deduction of the daily solar radiation 

in a tilted panel Eq.1; vi. Estimation of the maximum real water 
requirements Vmax using Eq. 15; 

iii. Calculation of the cell temperature 


which depends on characteristics of the vii. Estimation of the maximum flow 
panel Eq.4; required during the critical periode using 
Eq. 16; 
iv. Deduction of the panel yield using 
Eq.3; viii. Calculation of the pumping duration 
using Eq. 17 ; 
Gs npv Qmax At 


’ ' 1 


ix. Calculation of the PV panel surface Si using Eq. 18; 


x. Calculation of Ppv corresponding to the minimum panel surface using Eq. 19; 


| 
$ 


Sopt = sup(Si) 


t 


| Stop 


Figure 3. Main flowchart of optimal sizing of the photovoltaic pumping system for irrigation 


Data: a, 3, ET, and P, for the chosen site; 7),,8,,, and NOCT for the chosen panel; A’, for the 
chosen crop; 

Result: Optimal PV panel surface Sopt; 

for i= 1] : Nr do 


Estimation of Gs (&, 8) (i), Te(i) and ET,.(i); 
Deduction of Nps (i); 
if Pe(i) — ET.(i) < 0 then 
calculation of the total water requirements V,;(i) : Vii(i) -r : 
estimation of the maximum real water requirements volume: j 
Vinax (i) esum(V,; (i)) ro Vstorage (15) 
Vinax | 
estimation of Qmaz, the maximum flow required: Qmax iradi (16) 
min(t; 
P 
calculation of At, the pumping duration [1] :At =o (17) 
else 
pump water to storage tank: Vstorage <Nstorage(Vi + Vri); 
end 
At. mar 
calculation of the panel surface :S(7) ALO mas _ (18) 
Gr mp Me Npv 
end 
calculation of the photovoltaic power : Py, —SiG's(a, 3) Npv (19) 


deduction of S,,1, the optimal panel surface : Sopt <-sup(Si) ; 


Algorithm 1. Optimal sizing algorithm of the photovoltaic pumping system 
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5. STUDY CASE 

For a case study of the proposed algorithm, a specific simulation of an area of 0.5 hectares in the Fez- 
Meknes region, Morocco (latitude: 33° 53’36 ”North, longitude: 5° 32’50 ”West) cultivated tomatoes focus 
on Figure 1. For this reason, the tilt angle must be fixed at 33° and facing to the south. To use a PV pumping 
system to irrigate the crops and water livestock, given the significant solar radiation during the growing season 
in the site, unlike a complicated diesel engine. In the Fez-Meknes region, in general, farmers use furrow or drip 
irrigation, which is a principal method for irrigating tomatoes, thanks to its economic advantages in terms of 
saving water and increasing production. 

We have chosen the RECOM RCM-345-6MA mono-crystalline photovoltaic module for the character- 
istics of the panel, which are presented in Table 1. For the data input Table 2 contains climate data, photovoltaic 
pumping system, well, and crops data. The decade average daily values of solar radiation Gr, insolation t;, 
air temperature T,, evapotranspiraton HT, and total effective precipitation P, were obtained from CROPWAT 
program. T, cell temperatures were recalculated according to (4). The sizing algorithm was implemented dur- 
ing the growing season from the first decade of April to the last decade of August. Table 3 shows the water 
requirements for livestock watering. The water requirements of these animals vary little from one season to 
another, increasing slightly in summer and decreasing slightly in winter. We can assume that the daily water 
needs will be the same. 

Table 4 represents the real evapotranspiration /'T,.(2), the total required water quantities for irrigation 
(Vri), the panel yield npu, and the pumping duration At. Table 5 describes the results of the proposed algorithm, 
the value of the surface of the panel S; as well as the photovoltaic energy P,,, which are calculated for the 
three scenarios: 

— Scenario 1: optimal sizing using the proposed algorithm. 


— Scenario 2: optimal sizing using the proposed algorithm without storage tank, which means direct irrigation 
from the pump. 


— Scenario 3: sizing without optimisation. 


Table 1. Characteristics of the panel 


PV module RCM-345-6MA PV module type Mono-crystalline 
Peak power(Pmaz) 345 Wp Efficiency 17.78 % 
Voc 46.6 V Vmpp 38.1 V 
Isc 9.48 A Impp 9.06 A 
Operating temperature —40 + 85°C Coeff./Pmpp —0.40° C 
Coeff./ Voc —0.32 ° C Coeff./I sc 0.048 %° C 
NOCT 45°C Dimensions 1956*992*40 mm 


* Electrical data of PV module measured at STC (G's = 1000W/ m2, 
Te = 25°C, Air mass (AM) = 1.5) 


According to Table 2, the most critical month is July, with precipitation of 2 mm and evapotran- 
spiration of 5.84 mm/day or 175.2 mm/month. This month requires more water for irrigation, about 204.7 
mm/m. Hence, the dimensioning of the July system components is selected as the system dimensioning. 
The total water requirement is illustrated in Figure 4 and also the culture evapotranspiration and the effective 
precipitation for the growing cycle of the crop. The first scenario presents the results of our proposed algorithm. 
We obtained the surface of the panel S; as well as the photovoltaic power Ppy using (18) and (19). Hence, we 
can deduce that the optimum surface of the PV panel that satisfies the requirements of consumers during the 
observed period is Sopi = sup(S;) = 33.1 m?. 

The results obtained are compared to the two scenarios (using the proposed algorithm without having 
the tank and the other dimensioning without optimization) to verify the accuracy of the proposed dimensioning 
algorithm. As we see, in the same critical month of July, the necessary surface of the modules obtained for the 
second scenario is 38 m?, and for the third scenario is 93 m?. By comparing these results with the ones obtained 
using the proposed algorithm, it can be seen that the optimal values of the surface panel, obtained by the novel 
proposed of optimization algorithm are lower than those obtained by the compared scenarios. The precision of 
the obtained results depends mainly on all relevant elements, from the photovoltaic pumping system, the well, 
the microclimate, the soil, the crops, the storage tank, to the irrigation system. To obtain the optimum sizing 
of the irrigation PV pumping system, the proposed algorithm performs well under the conditions of the subject 
and satisfies the need in the observed period. 


A novel algorithm for optimal sizing of stand-alone photovoltaic ... (Asmae Hafian) 


1840 g ISSN: 2088-8694 


Table 2. Data input for the Fez-Meknes region 


Time stage i (Months) Parameters 
Gr(Kwh/m?/day) ti(hour) Ta(°C') ETo(mm/m)  Pe(mm) 

January 8.9 4.3 9.65 37.8 66.8 
February 12 5.5 10.7 51.6 70.5 
March 15.2 6.0 12.9 74.1 62.9 
April 18.9 7.0 15 96.3 64.5 

May 22.7 8.7 17.7 126.3 40 
June 24.2 9.5, 21.65 150 11.8 

July 25:3 10.6 25.15 175.2 2 

August 23.6 10.1 25.55 169.2 2 
September 19.1 8.5 22.55 129.6 11.8 
October 14 6.7 18.8 87.3 44.3 
November 9.9 5.0 13.9 52.5 67.5 
December 8 4.0 10.5 38.7 70.5 


Table 3. Water requirements of cattle herd 


Animal type No. of animals Daily requirement (gallon/day) 
Horses 2 2*15=30 
Sheep 20 20*2=40 
Milk producing cows 5 5*15=75 
calves 2 3*5=15 
Total 160 = 0.60 (m? /day) 


Table 4. Climatic and irrigation parameters and panel efficiency 


Time stage i (months) Parameters 

ETe(mm/m) Vu(mm/m) —_ mpo(%) _ Ath) 
April 58.3 1.7 13.48 5.23 
May 108 68.1 13.34 11.37 
June 171.6 159.8 13.34 12.42 
July 206.8 204.7 12.95 13.86 
August 123 121.8 12.93 13.21 


Table 5. Algorithm results for different scenarios 


Scenario 1 Scenario 2 Scenario 3 
Months 5; (m?)_— Pay (Kwh) Si2z(m?) — Paya (Kwh) Sis (m?) Pog (Kwh) 
April 8.52 4.87 19.1 10.91 30.6 17.90 
May 16.95 7.796 25.7 14.96 51.0 30.18 
June 24.74 10.42 32.7 18.69 70.0 40.01 
July 33.10 12.49 37.8 20.42 92.80 53.04 
August 30.5 12.08 36.5 20.01 90 51.44 


250 
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Figure 4. Total water requirements (V,.;), culture evapotranspiration ET and effective precipitation Pe chart 
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6. CONCLUSION 

Solar energy is one of the highest solar insolation rates. In this sense, the use of photovoltaic water 
pumping for irrigation or drinking water supply to improve the socio-economic conditions of users, mainly 
in remote areas, becomes an attractive solution. For this reason, the development of optimal solutions taking 
into account the Energy-Water duo remains a major asset to have a considerable impact on the dimensions of 
sustainable development. 

This paper presents a systematic algorithm for sizing the improved PV pumping systems, considering 
all relevant elements, from the PV pumping system, wells, local climate, soil, to the irrigation system. As 
a result of this algorithm, the optimum size of the PV pumping system was presented, to meet the hydraulic 
power requirements in the best possible way during the entire observed period. The photovoltaic pumping 
system was designed to fulfill the irrigation water requirements of tomatoes cultivated near the Fez-Meknes 
region location and marked out by an irrigated area of 0.5 hectares and also the livestock watering. That work 
verifies that the electrical power of the photovoltaic generator is relatively lower than that obtained without 
optimization. 
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